Universality of scaling and multiscaling in turbulent symmetric binary fluids 
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We elucidate the universal scaling and multiscaling properties of the nonequilibrium steady states 
(NESS) in a driven symmetric binary fluid (SBF) mixture in its homogeneous miscible phase in three 
dimensions (3d). We show, for the first time, via Direct Numerical Simulations (DNS) that structure 
functions of the velocity and the concentration gradient exhibit multiscaling in 3d and extended self- 
similarity (ESS). We also find that, in contrast to the well-known passive scalar turbulence problem, 
structure functions of the concentration show simple scaling. We propose a new shell model for 
SBF turbulence which preserve all the invariances in the ideal limit of the SBF equations and which 
reduces to a well-known shell model for fluid turbulence in the zero concentration field limit. We 
show that the shell model has the same scaling properties as the 3d SBF equations. Our combined 
results from our DNS of the SBF equations and shell-model studies consistently bring out the 
multiscaling of the velocity and concentration gradient fields and simple scaling of the concentration 
field. 

PACS numbers: 47.27.eb, 47.27.ek, 47.27.Gs 



Q 

U 
d 



> 

oo 
m 

(N 



X 



INTRODUCTION 

The scaling properties of correlation functions near a 
critical point in equilibrium statistical mechanics have 
been well understood over the past few decades. How- 
ever, understanding similar power-law scaling behaviours 
in structure functions in a variety of turbulent flows re- 
mains an open problem in nonequilibrium statistical me- 
chanics 0. In recent years, significant progress has been 
made in the study of equal-time structure functions in 
the turbulence of fluids, magnetohydrodynamics (MHD) 
and, most notably, passive-scalars Q. By contrast, for 
symmetric binary fluid turbulence, statistical studies are 
still in its infancy. 

To appreciate the context and the necessity for a sys- 
tematic study of the scaling properties of equal-time 
structure functions in such symmetric binary fluid mix- 
tures, it is important to recall some lessons from standard 
equilibrium critical phenomena 0, 0]. For a d dimen- 
sional spin system near a critical point, the equal-time, 
correlation function g, between two spins separated by 
the vector r (r = |r|), and its spatial Fourier transform 
g have power-law scalings : 



g(r;t,h) 
g(k;t,h) 



G{rt v ,hjt^) 



i-2+n 



G(k/t v ,h/t*) 
k 2-r, 



(1) 



where, t = (\T — T c \)/T c , T and T c are the temperature 
and the critical temperature, respectively, h = H/kBT c , 
H is the external field, ks is the Boltzmann constant, k 
is the wavevector, k = |k|, v, A, rj are critical exponents, 
and G and G are scaling functions. We note, in passing, 
that away from the critical point the correlation functions 



decay exponentially and the associated correlation length 
£ c diverges near a critical point as £ c ~ t~ v , if h = 0. 

Can we generalise such ideas of equilibrium statisti- 
cal mechanics to the case of homogeneous, isotropic tur- 
bulence in various settings? Indeed, the power-law be- 
haviours of equal-time structure functions, in the iner- 
tial range (to be defined later) , in fluid, passive-scalar or 
MHD turbulence have a certain similarity to the alge- 
braic dependence on r of correlation functions in critical 
theory. To make this connection explicit and lay the 
ground for our subsequent discussions, we begin with the 
increments of the longitudinal component of the velocity 
6u\\ (x, r, t) = [u(x + r, t) — u(x, t)} ■ (r/r), where u(x, t) 
is the velocity of the fluid at the point x and time t, 
and the subscript || implies the longitudinal component. 
The order-p, equal-time structure functions for the fluid 
(superscript u) field are defined conventionally as 



S«(r) = ([<Wx,r,i)n~rC : 



(2) 



the angular brackets indicate averages over the steady 
state for statistically steady turbulence or over statisti- 
cally independent initial configurations for decaying tur- 
bulence. The power law behaviour of such structure func- 
tions, which is valid for separations r in the inertial range 
r)d <C r -C L, where r\d is the Kolmogorov dissipation 
scale and L the large length scale at which energy is in- 
jected into the system are characterised by the equal-time 
exponents (£. 

Kolmogorov's phenomenological theory [l], [H, Q of 
1941 (K41) for fluid turbulence, predicts simple scaling 
(p' Kil — p/3- Subsequent experimental and numerical 
studies, however, strongly suggests the existence of equal- 
time multiscaling : Q is a nonlinear, convex, monotone- 
increasing functions of p. Indeed, it is important to 
remember that for the simplified stochastic Kraichnan 
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model [1, 0-Q of passive-scalar turbulence, multiscaling 
of equal-time structure functions can be demonstrated 
analytically. The analogue of the K41 theory for passive- 
scalar turbulence is due to Obukhov and Corrsin fiol fTll j . 
For the Schmidt number Sc = v/n ~ 1, where v is the 
kinematic viscosity of the fluid and n is the diffusivity of 
the passive scalar, the Obukhov- Corrsin theory leads to 
K41 scaling exponents for the passive-scalar case. 

In sharp contrast to fluid and passive-scalar turbu- 
lence, a systematic theoretical and numerical study of the 
statistical properties of symmetric, binary fluid (SBF) 
mixtures in three dimensions (3d) is still in its early 
stages and experiments performed on such systems have 
been typically concerned with measurements of effective 
transport coefficients [12j. Our prime concern here is to 
extend the ideas of equal-time scaling and multiscaling 
to the turbulence of SBF. In this paper we provide for 
the first time, via detailed Direct Numerical Simulations 
(DNS) and a new shell model that we propose for such 
a system, a systematic study of the statistical properties 
of equal-time, two-point structure functions in a statisti- 
cally steady, turbulent SBF mixture. We thus consider 
an incompressible, binary fluid mixture, with compo- 
nents labelled A and B, and having densities pa(*-, t) and 
/9b (x, t), respectively, such that the concentration field 
t) is defined via ip(x.,t) = |/»a(x, t) — ps(x, t)]/po, 
where po is the mean density. Furthermore, since we will 
be interested in a symmetric, binary fluid mixture we 
impose the constraint (ip(x.,t)) — 0. We elucidate the 
universal properties of homogeneous, isotropic SBF tur- 
bulence, in the absence of any macroscopic (mean) con- 
centration gradient, by measuring the scaling exponents 
of the equal-time structure functions of the velocity field 
u, the concentration field ip and the concentration gra- 
dient field b = Vtp. We show for the first time that al- 
though the exponents associated with b are multiscaling 
(like the exponents for u), the equal-time exponents for 
ip show simple-scaling. Our results are similar to the nu- 
merical quasi-Lagrangian (in two-dimensional flows) (l3j 
and in agreement with the predictions of one-loop field 
theoretical [H| studies of the SBF system. 



MODEL EQUATIONS 

In order to describe the coupled dynamical evolution of 
the field u and ip we need coupled dynamical equations 
for u and ip. The equation of motion of the velocity 
u is the generalised Navier-Stokes equation which now 
includes the stresses from the b field |15l Il6| 



where Ai, A2 > are coupling constants, and advection 
diffusion equation for ip: 



dip 

-jg + A 3 u • 



rjV 2 ip 



■ ft/>- 



(4) 



In Eqs.[3]and[H P and po are the local (effective) pressure 
and density, respectively; since we consider an incom- 
pressible fluid we further have po = const, and V • u = 0. 
The constants v and rj are the kinematic viscosity and 
concentration diffusivity, respectively. The functions f 
and /.0 are forcing terms which drive the system to a sta- 
tistically steady state. Galilean invariance of the system 
enforces Ai = A3 = 1 14, ltjj . Further, A2 may be set 



to unity by appropriately choosing the unit of ip (equiv- 
alently, by exploiting the rescaling invariance of ^)[14j|. 
Thus, in what follows, we set Ai = A2 = A3 = 1. It is 
clear from Eqs. © and Q (see also Refs. jlH. HE IbE]) 
that in the dynamics of a symmetric binary fluid mix- 
ture, the velocity field u couples with the concentration 
gradient b = and not with tp itself. Thus it is useful 
to write the coupled evolution equations in terms of u 
and b: The resulting equations are 

rhi VP 

— + (u-V)u = bVb + z/V 2 u + f (5) 

at po 

and the advection-diffusion equation for b is [l5|, [l6| 
dh 



dt 



V • (u • b) = -qV 2 b + g. 



(6) 



du 



+ Ai(u- V)u : 



VP 

Po 



■\ 2 VipV 2 ip + i>V 2 u + f, (3) 



Here, g = V/^,. Note that V x b = 0, thus b is an 
irrotational field. In a symmetric binary mixture, ip is not 
advected passively by the velocity field, but is active, i.e., 
the concentration gradient b reacts back on u and thus 
modifies the flow. Furthermore, since we are interested 
in the isotropic and homogeneous case, i.e., we have no 
mean concentration gradient, we impose (b) = 0. 



STRUCTURE FUNCTIONS AND 
MULTISCALING 

The order-p, equal-time structure function is defined 
as Sp(r) — (|a(x + r) — a(x)| p ), where a can be u, b or 
ip, x, r are spatial coordinates and the angular brack- 
ets represent an average over the NESS. For r in the 
inertial range which lies between the large length scale 
L and rjd, the Kolmogorov scale where dissipation be- 
comes significant, and at high fluid and concentration- 
gradient Reynolds numbers, Re and Reb, respectively, 
we expect power-law scaling 5f(r) ~ . The determi- 
nation of the exponents Cp has been one of the central, 
but still elusive, goals of studies in the statistical theory 
of turbulence. The extension of Kolmogorov's 1941 the- 
ory [B| to homogeneous, isotropic SBF turbulence, with 
no mean concentration gradient, yields Cp = p/3, i-e., 



3 



simple scaling. In isotropic and homogeneous pure fluid 
turbulence, we have corrections to simple-scaling expo- 
nents such that the equal-time exponents for such sys- 
tems (p = p/3 — 5Q, where SQ > and Cp is a nonlin- 
ear, monotonically increasing functions of p. Extensive 
analytical and numerical studies on the well-known pas- 
sive scalar problem [17| , which is the passive limit of the 
system considered here, clearly demonstrate that Cp has 
multiscaling qualitatively similar to Cp- I* 1 contrast, it 
has been shown for two-dimensional flows in Ref. fl3| . by 
using a Lagrangian approach, that Of = p/3, i.e., S^(r) 
shows only simple scaling. Ref. [14{ used symmetry ar- 
guments to show that Cp = p/3 and suggested that Cp 
should show multiscaling akin to Cp ■ It is thus expected 
that the multiscaling behaviour of the NESS in homoge- 
neous and isotropic SBF turbulence is characterised by 
Oj and Cp- Here, we confirm this in numerical studies of 
3dSBF equations and our shell model equations. 

Before we embark upon a discussion of our results, 
we explore the formal similarities between the dynami- 
cal equations of binary fluid turbulence and MHD. These 
become apparent when Eqs. ^ and (UJ) are compared 
with the incompressible 3dMHD equations. The incom- 
pressible 3e?MHD equations are given by [3] 

du . „, VP (V x B) x B „ 

— + (u ■ V)u = + + vV 2 u + f ; 

dt p a inpo 



A^k 



d-y\ 



<m 
~bt 



V x (u x B) = /i V 2 B + g. 



(7) 



Here B and /io are the magnetic field and magnetic 
viscosity, respectively. Other symbols have the same 
meaning as in Eqs. ([5]) and ©. The similarities be- 
tween Eqs. ©, © and |J7J| are noteworthy: (i) The 
concentration gradient field b and the magnetic field B 
have the same naive dimensions (which in turn is the 
same as the naive dimensions of the velocity u); and 
(ii) the non-linear terms in Eqs. ((5|) and (j6]) have ana- 
logues (in the sense of the number of field and gradients) 
in Eqs. (J7]) with same naive dimensions. All these sug- 
gest that the concentration gradient field b plays the 
role of the magnetic field in MHD 19]. In homoge- 
neous and isotropic 3dMHD turbulence, structure func- 
tions Sp{r) = (|a(x + r) — a(x)| p ) where a refers to both 
u and £?, exhibit multiscaling similar to pure fluid turbu- 
lence. This suggests that the structure functions S^.(r), 
like the magnetic field structure functions in MHD [20|, 
should exhibit multiscaling akin to pure fluid turbulence. 
In this paper, we confirm this conjecture. 

A promising starting point for a systematic the- 
ory is one where Eqs [5] and [5] are forced by Gaus- 
sian random forces f and g [cf. Refs. [2l[ for 
an application of this approach in pure fluid tur- 
bulence], whose spatial Fourier transforms, f(k, t) 
and g(k, t), respectively, have zero mean and co- 
variances (/i(k, t)fj(-k, 0)) 

(si(M)ff i (-k > o)) = A^h 



to noise variance (/^(k, t)f- l p(— k, 0)} 
where k is a wavenumber, t time, i,j Cartesian compo- 
nents in d-dimensions, Af and are a constant am- 
plitude and fy(k) = #y — kikj/k 2 is the transverse 
jrojector which enforces the incompressibility condition 
14j . One-loop renormalisation group studies 1J, [l6| of 
this model yield K41 energy spectra for u and b fields: 
E u ' b (k) ~ k 2 S^ b (k) ~ £r 5 / 3 for d = 3 and y = 4. Nev- 
ertheless, these RG studies have been criticised for a va- 
riety of reasons [22| such as using a large value for y in a 
small-?/ expansion and neglecting an infinity of marginal 
operators (if y = 4). These criticisms of the approxi- 
mations, however justified they may be, cannot be used 
to rule out the randomly forced model as an appropriate 
theory for SBF turbulence. In this paper, we use random 
Gaussian forcing in all our numerical simulations. 

We end this Section by summarising the key results 
in this paper. Our studies yield many interesting re- 
sults: The multiscaling exponents for u and b fields 
which we obtain from 3<iSBF and our shell models 
agree [Figs. [21 d] and Cp ~ Cp He close to, but be- 
low, the She-Leveque prediction (SL) Q for pure flu- 
ids (Cp 51, = p/9 + 2[1 - (2/3) p / 3 ]). Furthermore, the 
probability distribution functions (PDF) (Fig. [5]) for 
Sa a (r) = a Q (x + r) — a Q (x), a = u,b show non-Gaussian 
tails, whereas the same for 6ip(r) shows good agreement 
with a Gaussian distribution. These features of the PDF 
confirm the multiscaling behaviour of u, b and the simple 
scaling of tp. Earlier studies of fluid [24| and MHD [20| 
turbulence show that an extended inertial range is ob- 
tained if we use Extended Self Similarity (ESS): Thus, 
by making use of ESS, in which Cp/Cf follows from 
Sp ~ [iSg]^^ 3 , a — u,b, and ip, we expect, by anal- 
ogy, that the scaling range r extends down to r ~ 5rjd- 
We confirm this in our simulations for SBF turbulence. 



AfP^k^-ysit), 

v 5(t) (corresponding 



SHELL MODEL FOR THE SBF MIXTURE 

As is well known in turbulence, it is important to re- 
solve the large ranges of both temporal and spatial scales 
well. A DNS approach to hydrodynamical partial differ- 
ential equations, such as the one we use for the SBF mix- 
ture, is often very difficult if we want to resolve all the 
scales relevant to turbulence. To gain insight it is thus 
useful to consider simplified models of turbulence that 
are numerically more tractable than the partial differ- 
ential equations themselves. Shell models are important 
examples of such simplified models; they have proved par- 
ticularly useful for testing ideas of multiscaling in fluid, 
passive-scalar and MHD turbulence QMESIii. Keep- 
ing this in mind, we derive below a new shell model for 
the gradient of the concentration field in the SBF sys- 
tem and solve it numerically to obtain results in support 
of our DNS results. We should point out that for SBF 
turbulence, a shell model was derived in Ref. [27| for the 
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scalar concentration field and its coupling with the fluid 
field. 

Shell models cannot be derived from the hydrodynam- 
ical equations in any rigorous way. Such models are con- 
structed on a basis of a discretised Fourier space with 
logarithmically spaced wave vectors k n — k X n ,X > 1 
which are associated with shells n and dynamical com- 
plex, scalar dynamical variables which mimic, e.g., ve- 
locity increments u n over scales oc l/k n . Furthermore, 
we impose that k n be a scalar because spherical sym- 
metry is implicit in Gledzer-Ohkitani-Yamada (GOY)- 
type shell models which study homogeneous, isotropic 
turbulence 



261 ] . The logarithmic discretisation of the 
Fourier space allows us to reach very high Reynolds num- 
ber (which are impossible using DNS in present-day com- 
puters) even with moderate values of AT, where N is the 
total number of shells. 

The temporal evolution of such a shell model is gov- 
erned by a set of ordinary differential equations that have 
certain features in common with the Fourier-space ver- 
sion of the hydrodynamical equation. Thus, the shell 
model analogue of the Navier-Stokes equation, for ex- 
ample, will have a viscous-dissipation term of the form 
—vk\u n and nonlinear terms of the form ik n u n u n ' that 
couple velocities in different shells. (We note in passing 
that gradients appear as products of k n in shell models.) 
In the Navier-Stokes equation all Fourier modes of the ve- 
locity are coupled to each other directly but in most shell 
models nonlinear interactions are limited to shell veloci- 
ties in nearest- and next-nearest-neighbour shells. Hence 
sweeping effects common to equations of hydrodynamics, 
are absent in shell models. 

Keeping in mind the constraints of the hydrodynamical 
equations themselves, we propose the evolution equations 
for the shell model analogues of the velocity u n , the con- 
centration field ipn (see also Ref. 27 1), and the gradient 
of the concentration field b n as 



+ vkn]u n = i[Aik n u n+ iu n+2 + A 2 k n - 1 u n ^iu n+ i 



+ A 3 k n ^ 2 u n -iu n -2 + Aik n b n+ ib n+ 2 
+ A 5 k n b n+ ib n -i + A 6 k n2 b n -ib n - 2 }* 

+ fn, 




FIG. 1. (Color online) Log-log plot of the fluid kinetic energy 
E u (k n ) versus the wavenumber k n (blue filled circle joined by 
the continuous line) from our shell model studies. The thick 
black line indicates the K41 scaling. 



and 



l dt 



(8) 



+ V k l]K = i[A 7 k n (u n+ ib n+2 + u n+2 b n+1 ) 
+ A 8 fe„_i(u n+ i6 n _i + u„-i6„ + i 

+ A g k n - 2 {Un-lb n -2 + ti„_2&n-l)]* 

+ 9n, (10) 

respectively. In these equations, complex conjugation is 
denoted by *, and the coefficients are chosen such that 
the shell model analogues of total energy and the total 
autocorrelation of the concentration field is conserved in 
the absence of forcing and dissipation. Thus we obtain 
Ai = 1, A% = e - 1, A3 = e, ^4 + ^8 + ^9 = 0, A 5 - 
A 7 + A g = 0, A 6 + A 7 + A s = 0, A 7 + A g /X 4 = 0, 
A 7 - A 8 /X 2 = 0, and A$ + A 9 /X 2 = 0. We use the usual 
GOY model choice 25] of e = 0.5 and fix A 7 = 1 in order 
to obtain the values of the remaining constants. We have 
checked that our results are insensitive to the choice of 
A 7 . 



RESULTS FROM SHELL MODEL AND DNS 
STUDIES 



We begin by describing results from our numerical sim- 
ulations of the shell model for the SBF mixture. In our 

r "> 1 ; 2i ; -u 1 1 1 1 \ simulations the shell number is chosen such that 1 < n < 

[— + rjk n \yj n = i[k n (ip n +iu n -i — ipn — lu n +i) — — 

at N, where N = 22 is the total number of shells and we 

_ k n -i . ^ ^ ^ . use the boundary conditions u n — ijj n = b n = OVn < 1 

2 " " n n or Vn > Y. We use a second-order Adams-Bashforth 

k n +i t , , xl * method to solve the equations, a time step St = 10~ 4 

7T~ \Vn+2U n +l + tPn+lU n+2 )\ , , „_ 8 , ,. w , 

2 and v = rj = 10 m all our simulations. We choose a 

+ g n , (9) Gaussian, stochastic forcing on the fourth shell (n = 4) 
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FIG. 2. (Color online) Plots of with error bars from our 
128 3 DNS (blue triangle) and shell model (red star), C,f with 
error bars from our 128 3 DNS (blue diamond) and shell model 
(red filled-circle), and K41 scaling (thick black line) versus p. 
The lines connecting the data points from our simulations 
are a guide to the eye. The data from our DNS (upto 1 < 
p < 6) and shell model (shown for 1 < p < 6) are almost 
indistinguishable from each other upto p = 6. 
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FIG. 3. (Color online) Plots of Qj with error bars from our 
128 3 DNS (blue triangle) and shell model (red star) , and K41 
scaling (thick black line) versus p. The lines connecting the 
data points from our simulations are a guide to the eye. The 
data from our DNS (upto 1 < p < 6) and shell model (shown 
for 1 < p < 6) are almost indistinguishable from each other 
upto p = 6. 



to drive the system to a statistically steady state. Al- 
though in most studies of shell models, a deterministic 
force is used, we chose a stochastic forcing to make our 
shell model simulations consistent with our DNS. A snap- 
shot of the fluid kinetic energy spectrum, obtained from 
our shell model studies, which gives a good indication of 
the extent of the inertial range obtained, are shown in 
Fig. [I] with the K41 scaling indicated by the thick black 
line. The extent of scaling, a little over 3 decades, is typ- 
ical of such shell models and which allows measurements 
of scaling exponents with a higher degree of precision 



and confidence than in most DNS [261] . We show in Ta- 
ble 1 our equal-time scaling exponents if) (column 2), b 
(column 3) and u (column 4) fields; these exponents are 
calculated by using ESS, with respect to the third-order 
structure functions, for 50 different statistically indepen- 
dent statistically steady state configurations and quote 
the mean of these as our exponents and the standard 
deviation as the error-bars. We show the exponents for 
the velocity (red star) and the concentration field (red 
filled-circle), in Fig. [2] and for the gradient of the concen- 
tration field (red star), in Fig. G2 as a function of p; it 
is clear from the figures that there is clear multiscaling 
of the exponents associated with u (see Fig. [5] , red star) 
and b (see Fig. [3] , red star) and that the two agree with 
each other within error-bars (compare columns 3 and 4 
in Table 1). In contrast, the exponents for ip (see Fig. [2] 
, red filled-circle) shows simple scaling and is indistin- 
guishable, within error-bars, from the K41 prediction. 



How much of our results from the shell model studies, 
as described above, carries over to the actual Direct Nu- 
merical Simulations of the SBF mixture? We now de- 
scribe the results obtained from our pseudospectral stud- 
ies of the randomly forced symmetric binary fluid mix- 
ture equations (Eqs[3] and 2]) and compare them with 
those obtained from the numerical solutions of our shell 
models. In our DNS, we keep y = 4, corresponding to 
K41 spectra for the u and b fields, we use resolutions 
of 96 3 and 128 3 with a cubic box of linear size L = 2tt, 
and periodic boundary conditions. Our numerical scheme 
is identical to that in Ref. 28]. We use hyperviscos- 



ity and hyperdiffusivity together with ordinary viscosity 
and diffusivity. For the resolution 128 3 , we are able to 
achieve Taylor microscale Reynolds number Re\ ~ 150. 
In Fig HI we show a log-log plot of the fluid kinetic en- 
ergy spectrum obtained from our DNS. Although our 
Reynolds number is not very high, we do obtain an in- 
ertial range close to three-quarters of a decade as can 
be seen in the figure. In the NESS obtained from these 
DNS we calculate the exponents by using ESS with 
respect to the third-order structure function, from log- 
log plots of Sp(r) versus <Sf(r) (a = u i7 bi 7 ip). From 
such plots, we use a modified local slope approach to 
obtain the equal-time exponents : We calculate the ex- 
ponents over various ranges within the inertial range; we 
quote the mean as our exponent and the standard devi- 
ation as the error-bar. We find that : (i) The exponents 
Cp (Fig. H blue triangle) and (p (Fig. ED blue triangle) 
display multiscaling very similar to that in fluid turbu- 
fence:' C?/(? > 2/3, Cp/CT < Pl^P > 3, m = u, 6, 
and, Q/Q an d Cp/Cl are equal to each other within our 
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error-bars (compare columns 6 and 7 in Table 1); and (ii) 
Cjf/Cl ~ p/3 (Fig.H blue diamond). In addition, we cal- 
culate the normalized probability distribution function 
(PDF) P[Sa(r)} (a = u^b^ip) for r/rj d = 7.7 (see Fig. 
[5]). We find P[8u(r)] and P[Sb(r)] are nearly overlapping 
and have much longer tails than P[8%p{r)\. Furthermore, 
P[S4>(r)} is well represented by a Gaussian of unit vari- 
ance. In Fig. [6] we show plots of P[(5 - !/ , ( r )] versus 8ip(r) 
for three different separations r/r/d = 7.7, 23.2, 26.2 in 
the inertial range. A Gaussian of unit variance is again 
shown for comparison. We find that for all values of r, 
the plots overlap with each other and with the Gaussian. 
Also similar PDF plots for u (Fig. [7] and b (Fig. [5J for 
three different separations r = jr\d = 7.7, 23.2, 26.2 not 
only show a marked departure from a Gaussian (as in- 
dicated by a continuous dark blue line) as was seen in 
Fig. [5] but also no collapse of the curves for different r 
(unlike the case for ip). These PDFs further strengthen 
and provide compelling evidence for our main results (i) 
and (ii) above. We present, in Table 1, the multiscaling 
exponents <Qf (column 5), (column 6), and ip™ (col- 
umn 7) from our DNS. We have checked that our results 
from the two different resolutions for the DNS agree with 
each other within error bars. 

Given our modest resolution for our DNS, it is useful 
to examine how far we are justified in calculating mo- 
ments up to order 6. It is well known that for higher or- 
der moments, the large contributions from the tails of the 
PDFs make statistical convergence progressively poor. In 
order to study statistical convergence, a good prescrip- 
tion is to examine the convergence of the moments of the 
differences of the velocity and the concentration fields 
29]. Hence we study the bulk contributions Cq [8a] to 
the sixth-order structure function (the highest order for 
which we present results in this paper) Sg(r) (a = ip, u) 
defined as: Cq(x) — J Q x 6 p(x) where x = 8ip(r) or 8u (all 
suitably normalised), and p(x) is the PDF of the same, 
respectively. In Fig. HI we show a plot of Ce[8ip(r)] ver- 
sus 8ip(r) for two different r = 26. 2r)d (black) and 7.7^ 
(red): The two curves overlap, as is expected for a Gaus- 
sian form for P[Sip(r)] for various different values of r. 
In Fig. [10] we show a similar plot for C§[8u{r)\ for the 
same two r as before. Due to the non-Gaussian nature of 
P[8u{r)\, the two curves for two different r do not over- 
lap. These plots strongly display statistical convergence 
of the corresponding sixth-order moments. We obtain 
similar convergence for the gradient of the concentration 
field which we do not show here. 




FIG. 4. (Color online) Log- log plot of E u (k) versus k/k d (blue 
*) from our DNS studies (see text), kd ~ 98 is the dissipation 
scale wavenumber here. The thick black line indicates the 
K41 scaling. 




B| 



FIG. 5. (Color online) Semilog plots of the probability distri- 
butions P[Su(r)], P[8b(r)] and P[5tp(r)] versus r in the iner- 
tial range, averaged over both time and the Cartesian compo- 
nents, from our DNS; a Gaussian distribution (blue continu- 
ous line) is shown for comparison. 



CONCLUSION 

In summary, then, in this paper we have investigated 
the scaling and multiscaling properties of a turbulent 
symmetric binary fluid mixture via detailed numerical 
simulations. We find that Sp ,b (r) exhibit multiscaling 
similar to fluid turbulence in the inertial range, whereas 
S^(r) exhibit simple p/3 K41-scaling (within error bars). 
Moreover, the probability distributions P[8u(r)] and 
P[8b(r)] are nearly overlapping and have tails longer than 
that of P[8tp(r)] for r in the inertial range. We also pro- 
pose a new shell model for the gradient of the concen- 
tration field and numerically solve it as well as the shell 
model for scalar concentration field. Our results from 
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FIG. 6. (Color online) Semilog plots of normalized proba- 
bility distributions P[5tp(r)] as functions of 5ip(r) for three 
different separations r in the inertial range. A semilog plot of 
normalized Gaussian is shown for comparison. All plots are 
overlapping (see text). 








FIG. 8. (Color online) Semilog plots of normalized probability 
distributions P[Sb(r)] as functions of Sb(r) for three different 
separations r in the inertial range. A semilog plot of normal- 
ized Gaussian is shown for comparison (see text). The values 
of r increases going inwards. 




FIG. 7. (Color online) Semilog plots of normalized probability 
distributions P[8u(r)\ as functions of 8u(r) for three different 
separations r in the inertial range. A semilog plot of normal- 
ized Gaussian is shown for comparison (see text). The values 
of r increases going inwards. 



by the convincing PDFs that we obtain for various quan- 
tities in our DNS and the clear evidence of statistical 
convergence which justifies measurements of equal-time 
exponents upto order 6 in our DNS. 

Our results may be explained from the analytical 
framework based on symmetry arguments developed in 
Ref. [14j, where it has been shown that the presence 
of an additional continuous symmetry (kind of a gauge 
symmetry), not present in the passive scalar turbulence 
model, is responsible for the simple scaling behaviour of 
Cjf . It would be interesting to investigate the properties 
of the turbulent NESS of SBF at low temperature, below 
the consolute point, when instabilities leading to phase 
separation competes with turbulent mixing. Work is in 
progress in this direction. Finally, our results may be 
tested in experiments similar to Ref. |3Q |. One of the 
authors (AB) gratefully acknowledges MPG(Germany)- 
DST(India) for partial financial support through the 
Partner Group program. 



our shell models are in agreement with our DNS stud- 
ies. Our results are 3d analogues of those of Ref. [l3j ]. 
where simulations with particles in 2d flows were used 
to show that the structure functions of the concentration 
field for the SBF problem do not multiscale. The re- 
sults from both our DNS of the SBF equations and shell- 
model studies are complementary to and agree well with 
each other and bring out the multiscaling of the velocity 
and concentration gradient fields and simple scaling of 
the concentration field. The validity of these conclusions 
is strengthened not only by the reliability of the scal- 
ing ranges usually associated with the measurement of 
equal-time structure functions in shell models, but also 
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